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Abstract 

We study the motion of gap solitons in two models of apodized nonlinear fiber Bragg gratings 
(BGs), with the local reflectivity k varying along the fiber. A single step of k, and a periodic 
array of alternating steps with opposite signs (a "Bragg superstructure") are considered. These 
structures may be used in the design of various optical elements employing the gap solitons. A 
challenging possibility is to slow down and eventually halt the soliton by passing it through the step 
of increasing reflectivity, thus capturing a pulse of standing light. First, we develop an analytical 
approach, assuming adiabatic evolution of the soliton, and making use of the energy conservation 
and balance equation for the momentum. Comparison with simulations shows that the analytical 
approximation is quite accurate, unless the inhomogeneity is too narrow, or the step is too high: 
the soliton is either transmitted across the step or bounces back from it. If the step is narrow, 
systematic simulations demontrate that the soliton splits into transmitted and reflected pulses 
(splitting of a BG soliton which hits a chirped grating was observed in experiments). Moving 
through the periodic "superstructure" , the soliton accummulates distortion and suffers radiation 
loss if the structure is composed of narrow steps. The soliton moves without any loss or irrieversible 
deformation through the array of sufficiently broad steps. 

PACS numbers: 42. 81. Dp; 42.50.Md; 42.65.Tg; 05.45.Yv 
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I. INTRODUCTION 



Bragg gratings (BGs) are periodic structures produced by a periodic variation of the 
refractive index along a fiber or an optical waveguide. Devices based on fiber gratings are 
widely used in optical systems Gap solitons (in this context, they are also called BG 
solitons c 2jj]) are supported by fiber gratings through the balance between the BG-induced 
linear dispersion (which includes a gap in the system's linear spectrum) and Kerr nonlinearity 
of the fiber material. Analytical solutions for BG solitons in the standard fiber-grating model 
are well known P, 0] . Studies of the stability of these solutions have led to a conclusion that 
almost exactly half of them are stable (see details below) M. These solitons were created 

HQ 

and studied in detail in the experiment |a UJ . 

Recently, much attention has been attracted to the creation of ''slow light" in various 
media including fibers (for instance, a thin fiber embedded in an electromagnetically 
controlled molecular solid 9]). In particular, the possibilities to generate slowly moving 
optical solitons were considered in various settings The nonlinear fiber grating is a 

medium where it may be possible to stop a soliton. Indeed, solutions for zero-velocity gap 
solitons are available, in which the counter-propagating waves are mutually locked in the 
dynamical equilibrium between the linear Bragg-resonant conversion and Kerr nonlinearity 
keeping them together through the XPM (cross-phase- modulation) interaction P, |^ . 

Actually, the BG solitons that were created in the experiments to date were fast ones, 
the slowest among them having the velocity equal to half of the light velocity c in the 
fiber [7]. A possibility to capture a standing BG soliton is to bind it on an attractive 
local inhomogeneity [rj . Q] (earlier, it was demonstrated that a local defect in the fiber 
grating can also stimulate a nonlinear four-wave interaction without formation of a soliton 



14j). Moreover, it is possible to combine the attractive inhomogeneity and local linear gain, 



which is a way to create a permanently pinned soliton in the presence of loss [15]. Stable 
pinned-solitons are not only objects of fundamental interest, but they also have an obvious 
potential for design of various nonlinear-optical elements. 

Another approach to the creation of standing BG solitons was explored in a recent work 
where it was demonstrated that head-on collisions between two moving solitons may 
lead to their fusion into an immobile soliton (with residual internal oscillations, which is 

n 

explained by the existence of an oscillatory intrinsic mode in the stable BG soliton |5j), 
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provided that their initial velocities are smaller than fa 0.2 c. Thus, one should still bridge 
the remaining gap between the experimentally available minimum soliton velocity, 0.5c, and 
the necessary value of 0.2c. 



A possibility to resolve this problem was also proposed and briefly considered in Ref. 11 1 
(a similar possibility was mentioned still earlier in Ref. passing the BG soliton through 

an apodized fiber grating, with the Bragg reflectivity gradually increasing along the fiber. In 
this case, the soliton may be adiabatically slowed down. An objective of the present work 
is to consider this problem in a systematic form, both analytically and numerically. Besides 
that, we also aim to study the motion of the soliton through a periodic array of alternating 
reflectivity steps with opposite signs. 

It should be mentioned that apodized fiber gratings (a gradually varying reflectivity may 
also be featured by chirped gratings [3, ill are widely used in various applications 
Some experimental [?], 13 an d theoretical |l7| results concerning passage of BG solitons 
through apodized gratings were reported earlier. These include the use of the apodization 
to facilitate the launch of solitons into the fiber grating observation of splitting of a BG 
soliton which hits a chirped grating Jlfil ]. and analysis of the soliton's dynamics in the case 
when the apodized-grating model reduces to a perturbed nonlinear Schrodinger equation 
jl7|. The formulation of problems considered in this work, and methods developed for the 
analysis, are essentially different from the earlier works, as will be seen below. 

The rest of the paper is organized as follows. The model is formulated in section 2, and 
in section 3 we work out the analytical approach, which is based on treating the gap soliton 
moving in the apodized BG as a quasi-particle obeying an equation of motion derived from 
the balance equation for the soliton's momentum. This adiabatic approach is relevant when 
the Bragg reflectivity varies on a spatial scale which exceeds the size of the soliton (for 
slow solitons, a physically relevant size is, typically, ~ 1 mm, while for fast solitons the 
size is additionally subject to the "Lorentz contraction", see below). Results of systematic 
numerical simulations of the full model with a step of the Bragg reflectivity (that will also 
be sometimes called "barrier" ) are presented in section 4. First, we compare the analytical 
predictions produced by the adiabatic approximation with numerical results, and conclude 
that the adiabatic approximation is quite accurate, provided that the step is not too steep. 
Further, we present numerical results for the case of steeper apodization. In that case, 
the soliton impinging upon the step may split into two secondary solitons (transmitted 
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and reflected ones), or get completely destroyed. In section 5 we briefly describe results of 
simulations of the motion of the soliton through a periodic array of alternating positive and 
negative steps. It is found that the soliton distorts itself and emits radiation if the steps are 
narrow. If they are broader, the deformation of the moving soliton accumulates slower, and 
it can move without any irreversible damage through the array of sufficiently broad steps. 
Section 6 concludes the paper. 



II. THE MODEL 

The commonly adopted model of the nonlinear fiber grating is based on a system of 
coupled equations for the local amplitudes of the right- (u) and left- (v) traveling waves , 



iu t + iu x + kv + (\u\ 2 /2 + |t>| 2 ) u = 0, 

iv t - iv x + ku+ (M 2 / 2 + M 2 ) v = 0, (1) 

where x and t are the coordinate and time, which are scaled so that the linear group velocity 
of light in the fiber is 1, and k is the Bragg-reflectivity coefficient. In the apodized grating, 
k is a function of x. 



Exact solutions to Eqs. (£[]) with k = const, which describe so. 



c (c 2 < 1) through the uniform BG, were found in Refs. j^] and 



itons moving at a velocity 



u sol = 1 /^^(l-c 2 ) 1/ V(X)exp [i<j>(X)-iT cos 9] 



v ^ = ~V 3-/ (l-c 2 ) 1/ V*(X)exp [i^-iT cosO]. (2) 

Here, the asterisk stands for the complex conjugation, 9 is an intrinsic parameter of the 
soliton family which takes values < 6 < it, and 

X = k (1 - c 2 )~ V2 [x - f (t)] , T = « (1 - c 2 y 1/2 (t - cx) , 

(j)(X) = -^-^tan" 1 {tanh [(sin 9)X]tan(9/2)} , (3) 
3 c 

W(X) = (sin 9) sech [(sin 9)X - % (9/2)} , 

where d£/dt = c (note the above-mentioned " Lorentz contraction" of the soliton's size, 
obvious in the definition of X). 



Equations conserve three dynamical invariants: the norm (frequently called energy in 
optics) and momentum, 



E= (\u\ 2 + \v\ 2 ) dx, P 



(u* x u + v* x v) dx, 



(4) 



and also the Hamiltonian (which will not be used below). For the soliton solution (J2J), the 
norm and momentum take values 

86(1 -c 2 ) 



E so \ = 



sol 



8kcVT 



% (sin 9 - 9 cos 9) + r 



(5) 



L(3-c 2 

(note that the energy does not depend on k). 

The expression for the momentum strongly simplifies in the "nonrelativistic" limit (c 2 <C 
1), which makes it possible to identify an effective soliton's mass, 



M sol = lim — = — (7sinfl-4flcosf 

c^o c 9 



(6) 



Note that the gap solitons are stable in the region 9 < 9 cr , where 6 CT is very close to tt/2 
[for instance, 9 cr ~ 1.01 (vr/2) if c = 0] p. In all the stability region, the mass (jHI) increases 
with 9, i.e., larger 9 corresponds to a "heavier" soliton. 



III. THE ADIABATIC APPROXIMATION 

In the apodized BG, with k = k(x), the momentum is no longer conserved; instead, the 
following exact equation can be derived from the underlying equations (JH): 

dP f +co c/k . , 

— = 2 / —Re (uv*) dx. (7) 

The adiabatic approximation applies to the case of a slowly changing k(x), when it varies on 
a length which is essentially larger than the soliton's size. As it follows from the expressions 
(J2J) for the soliton solution, this implies a general condition, sin# ^> y \ — c 2 (k _1 \dK/dx\). 
For slowly varying k(x), we rewrite Eq. ((Zj) as 

°±- = 2^ Re(uv*)dx, (8) 
dt dxj^ v ; ' y J 

assuming that the value of dn/dx is taken at the point x = £(£), where the soliton's center 

is located at a given moment of time, see Eqs. ©• Then, the left-hand side of Eq. (jHJ) 
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can be calculated with the unperturbed soliton solution (J2J), and the expression for the 
soliton's momentum from Eqs. (0) can be substituted in the left-hand side. This yields the 
following general evolution equation for the soliton's parameters [the velocity c(t) and the 
mass parameter 9(t)\: 



d 
dt 



;(OcVT 



(7 - c 2 ) 
(3 - c 2 ) 2 



(sin 9 — 9 cos i 



9 cos 9 
3-c 2 



dn (1 



c 2 ) 3/2 



sin 9. 



di 3-c 2 (9) 
Note that, as k(x) is taken at the point x = which changes in time according to 

d^/dt = c, the coefficient k is subject to the t-differentiation on the left-hand side of Eq. 

©■ 

Unlike the momentum, the energy of the wave fields remains the dynamical invariant in 
the presence of the apodization. Therefore, using the expression for the soliton's energy 
from Eqs. (|5|) . in the adiabatic approximation one can eliminate 9(t) in favor of c(t): 

3-c 2 E 
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(10) 



1-c 2 8 

Thus, Eqs. (JjJJ) and (JTUJ), together with the above-mentioned relations k = k(£) and d^/dt = 
c, constitute a closed system of the evolution equations for 9{t) and c(t). 

The adiabatic approximation strongly simplifies in the above-mentioned "nonrelativistic" 
limit, c 2 C 1: Eq. (|10|) then amounts to 9 = const, the remaining equation for c(t) being 

3 sin 9 dn 



It {KC) 



(11) 



7sin^-4^ cos9d£ 

Further, using the identities d(nc) /dt = [d(nc) /d£\ (d£/dt) = cd(K,c) jdx (as d^/dt = c 
we transform Eq. 



d (kg) 

which can be integrated to yield 

3 sin 9 



3 sin^ 



dn 



7sin^-4^ cos9d£ 



+ c 2 



3 sin 9 



k 2 (0 7sin^-4^cos^' 



(12) 



7sm9 - A9 cos9 

where the subscript refers to initial values. In particular, Eq. ()12j) implies that the soliton 
may be brought to a halt (c 2 = 0) at a point £hait where k 2 attains the value 

2 (7 sin 9 - 49 cos( 



K 



halt 



Kr 



3 sin^ 



(13) 
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Of course, the expression (fT3j) makes sense if ft^ alt does not exceed the largest value of k 2 
available in a given apodization profile. 

Note that, except for the special case when k 2 (£) = fi^ait is attained at £ = oo, Eq. (JT2J) 
predicts that the soliton will not be stuck forever at the halt point £ halt; but will actually 
bounce back. Indeed, with regard to the relation c = d^/dt, it follows from Eq. (|12j) that, 
around the halt point, the soliton's law of motion takes the form £hait — £ ~ (t — thait) 2 ? 
where thait is the moment of time at which the velocity c(t) vanishes. 

We also notice that, if the soliton passes an apodized region (step) across which k(x) 
increases, then the final value of the velocity, as given by Eq. (fT2|) . increases with 9. This 
complies with the above conclusion that the soliton is heavier for larger 9: due to the inertia, 
a heavier object suffers smaller deceleration passing a potential step. 

IV. NUMERICAL RESULTS 

The underlying equations ((TJ) were simulated with the step-wise profile of the local reflec- 
tivity, 

k(x) = + -A/ttanh ^- — — ^ , (14) 
where the step's height An is positive, and a normalization condition is imposed, 

K (= +oo) = k q + (Ak/2) = 1. (15) 

The width of the step is u>, and its center is located at the point The split-step 

fast-Fourier-transform method was used to run the simulations. 

In the simulations, the gap soliton was launched, as the exact solution (J2J) with a positive 
velocity Cq, from the left edge of the integration domain. Results will be presented, chiefly, 
for stable solitons, which implies that the initial values of 9 should be taken from the interval 
$o < tt/2 (see above). In fact, simulations were also performed for solitons with 9 > 7r/2. 
The results are not drastically different for them (a brief description of this case is also given 
below); in particular, the passage of the step does not catalyze the onset of the soliton's 
instability, which develops on a larger time scale. 
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A. Verification of the adiabatic approximation 



First of all, predictions produced by the above adiabatic approximation were checked 
for the case of a sufficiently smooth step (w not too small). The checkup was done in 
two stages. In the lowest approximation, the explicit solution (|12|) . that was derived in 
the "nonrelativistic" approximation, was used, along with the corresponding condition 9 = 
const. To obtain a more accurate result at the second stage, the full evolution equation Q, 
combined with the relation ()10|). was solved numerically. It was concluded that the adiabatic 
approximation is quite accurate, as long as the step remains smooth. This can be illustrated 
by typical results obtained for the soliton with the initial value of its intrinsic parameter 
$o = 7r /2 impinging on the step with the width w — 6. 

In this case, the special value of the initial velocity (co) halt , which provides for the halt 
of the soliton in the direct simulations, was sought for at fixed values of k,q (the "halt" was 
realized so that the soliton's velocity would drop to zero, and after being stuck for a long 
time, the soliton would then very slowly start to move backward). It was thus found that 
the numerical values are, for instance, 



On the other hand, for these values of cq the adiabatic approximation predicts which must 
be the final values of the reflectivity, K& na i = + (1/2) Ak [see Eq. (fill )], that will halt the 
soliton. In particular, in the "nonrelativistic" approximation, we simply set c 2 = in Eq. 
(JT2J), which yields 



The predicted value of Kfi na i must be compared with Kfi na i = 1, which was imposed by the 
normalization (|TK|) adopted in the full numerical simulations. 

The "nonrelativistic" approximation predicts, for the two cases indicated in Eq. (fT7)|) . 
ftfinai — 0.971 and 0.755, respectively. The former value is quite close to the exact one, 
ftfmai = 1- The latter value is not very close to it, but the corresponding initial velocity, 
Co = 0.74, does not correspond to the "nonrelativistic" case. On the other hand, the use 
of the full adiabatically derived equation © predicts, for the same two cases, na na \ = 1.017 
and 1.023, respectively. 

Systematic comparison of the numerical and analytical results is presented in Fig. CJ 



(co)hait («o = 0.8) = 0.45; (c ) halt (k = 0.5) = 0.74. 



(16) 




(17) 
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which shows the border between the reflection of the moving gap soliton from the apodization 
step, and its passage through the step, in the plane of (cq,A«). For smaller values of the step's 
size, the analytical approximation provides for results which are quite close to the numerical 
findings. Making the perturbation stronger (increasing An), we observe a growing deviation, 
but the full "relativistic" approximation, based on Eqs. (JOJ) and (jlUj) . still provides for a 
reasonable approximation even for larger An (and larger velocity), while the oversimplified 
"nonrelativistic" approximation, based on the single equation (j!2|) . becomes irrelevant in 
that case. 

B. Different outcomes of the collision of the gap soliton with the apodization step 

To present direct numerical results in a systematic form (especially, for the situation 
with a steeper apodization profile), we first fix the height and width of the step, setting 
Ak = 0.5 and w = 6 [note that this width is much smaller than it was in the case of Fig. Q 
(w = 20), hence the step is much steeper in the present case]. Launching the solitons with 
various initial values Co and 8q of the velocity and mass parameter to hit the apodization 
step, it was found that, as it might be naturally expected, more energetic solitons, with a 
large velocity c and/or larger mass parameter 6, pass the step with deceleration, a typical 
example of which is shown in Fig. 2(a), and slower (smaller c) and/or lighter (smaller 9) 
solitons bounce back from the step, see a typical example in Fig. 2(b). Note that the only 
difference between the cases displayed in Figs. 2(a) and 2(b) is that Co = 0.75 in the former 
case, and Cq = 0.7 in the latter case, showing that the border between the transmission 
and repulsion is in between these values. On the other hand, the calculation based on the 
quasi-particle equation of motion © predicts that the border is at c = Cq ~ 0.61 for the 
same case, which agrees reasonably well with the numerical results, considering that the 
apodization profile is rather steep, and the apodization step is high. The applicability of the 
adiabatic approximation to the present cases is also corroborated by the observation that 
very little radiation loss was generated in the direct simulations. 

If the bounce takes place at a point £hait located far to the right from the step, the soliton 
gets stuck there for a finite but long time, also in agreement with the prediction of the 
quasi-particle approximation. An example of such a case (which, obviously, is of special 
interest to the experiment and applications, suggesting a real possibility to capture a pulse 
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of "standing light") is shown in Fig. 2(c). Note that, in comparison with a more generic 
case of the ricochet shown in Fig. 2(b), in the case of the quasi-trapping of the gap soliton, 
the bounce indeed occurs after the soliton has advanced much farther to the right. 

In the same situation as considered above, but with smaller values of 9q, the quasi- 
particle approximation does not apply, as the soliton size becomes much larger than the 
step's width, hence the step may not be considered as a smooth one in any approximation. 
The simulations indeed produce a drastically different result in this case: hitting the step, 
the soliton splits into two pulses, transmitted and reflected ones, which is illustrated by Fig. 
2(d). In this case, both pulses eventually decay into radiation, rather than self-trapping into 
small-amplitude solitons. 

The controllable splitting of the incident soliton is an interesting effect in its own right, 
and it has a potential for various applications. Splitting of the BG soliton hitting a chirped 
fiber grating has already been observed experimentally and reproduced in numerical simu- 
lations Q|. In fact, in that case the incident soliton split into three pulses: a transmitted 
soliton and two reflected ones. 

In this connection, it is relevant to mention that, although we do not display detailed 
results for the (generally) unstable solitons with 9 > tt/2, the interaction with the step was 
simulated for them too. In most cases, the intrinsic instability of the soliton does not set it in 
the course of the limited time before it hits the step (which corresponds to the experimental 
situation, in which the segment of the fiber grating before the apodized region is not going 
to be very long [6]). Then, if the step is smooth, the heavy soliton passes it or bounces back, 
without catalyzing the onset of the intrinsic instability, and without conspicuous emission 
of radiation. The soliton emits an appreciable amount of radiation if the barrier is steeper. 
The splitting of the soliton with 9 > 71/ 2 into forward and backward moving parts, similar 
to what is shown for smaller 9 in Fig. 2(d), does not occur on the barrier with w = 6, which 
is the value selected for Fig. 2, suggesting that the soliton is a more cohesive object for 
larger 9 (see further discussion of this point below). Splitting takes place if the barrier is 
made still steeper and taller - for instance, with w = 4 and Ak = 0.7, see Fig. 3. As well 
as in the case displayed in Fig. 2(d), in the latter case the two pulses do not eventually 
self-trap into secondary solitons, but rather decay into radiation. 
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C. Scanning the parameter space 

Having displayed typical examples of different outcomes of the collision, we proceed to 
summarize the results in a systematic form. First, a comprehensive description of the tran- 
sition from the bounce of the soliton to the passage, together with the splitting (if any), are 
provided, for different values of the barrier's height Ak, by the plots in Fig. 4. For fixed 
9, they show the share e of the initial soliton's energy which is reflected back as a result of 
the interaction of the soliton with the step. A steep drop from e = 1 to e = in Figs. 4(b) 
and 4(c) implies the transition from the bounce to the transmission without splitting. On 
the contrary, a gentle crossover means that the impinging soliton was split into two parts, 
with the energy ratio between them depending on the initial velocity. These plots clearly 
show strengthening of the soliton's integrity (its stabilization against the splitting) with the 
increase of its energy. 

Further, in Fig. 5 we show the velocities Cf and q, of the transmitted ("final") and re- 
flected (bouncing) fragments of the incident soliton, together with the backscattered-energy 
share e, as functions of Ak. In the ranges where only one velocity is shown in Fig. 5, the 
soliton is transmitted or reflected as a whole, while the splitting takes place in intervals of 
the values of An where both continuous curves (for c/ and q>) coexist in Fig. 5. 

Note that the transition from the full transmission to full reflection, i.e., from e = 
to e = 1, takes place in exactly the same interval where the splitting occurs, cf. Fig. 4. 
As is seen, the splitting interval is very narrow for the heavier soliton, with 9 = 0.57T, and 
much broader for the lighter one, with 9 = 0.37T, which is another common feature with the 
dependences displayed in Fig. 4. 

The dependence of the outcome of the collision on the width of the apodization region 
was studied too, as it helps to understand how really wide the barrier must be to secure the 
adiabatic character of the interaction, and how narrow it is when the splitting occurs. In 
Fig. 6, we show the dependence of the final velocity of the transmitted soliton, c/, and the 
reflected-energy share, e, on the width while the step's height is fixed, Ak = 0.5. The final 
velocity of the transmitted soliton approaches an asymptotic value for large values of w. 

Figure 6 shows that the transition from the strongly non-adiabatic regime (with splitting) 
to a nearly adiabatic one is steep itself, taking place in a relatively narrow interval around 
w = 3. A practically important implication of the results presented in Fig. 6 is that, for the 
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experimentally relevant case, with the intrinsic size of the BG soliton ~ 1 mm the step 
which may be regarded as a sufficiently smooth one should have the width > 3 mm. 

V. DYNAMICS OF THE GAP SOLITON IN A PERIODICALLY- APODIZED 
STRUCTURE 

To test the robustness of the gap solitons in the fiber gratings with a variable Bragg 
reflectivity k, and its potential for the use in more sophisticated devices, we also simulated 
the motion of the soliton in the grating subjected to a periodic modulation of k. It was 
built as a periodic concatenation of steps with alternating signs of Ak; the fiber grating 
of this type may be considered as a Bragg "superstructure". This scheme can be realized 
either directly in a long periodically modulated grating, or as an apodized BG written in 
a fiber loop, although the stability problems in these two settings are not equivalent. The 
soliton dynamics in such a superstructure is a problem of interest in its own right - cf. the 
study of solitons in other periodic strongly inhomogeneous nonlinear optical media, such as 
ibers with dispersion management various schemes with "nonlinearity management" 
rg| (including layered media with a periodically changing sign of the Kerr nonlinearity 
"tandem" structures [2]]], the "split-step" model [22 1. waveguide-antiwaveguide concatena- 
tions , etc. A unifying feature of these systems is the surprising robustness of the solitons 
in them. 

Typical results produced by the simulations of the soliton propagation in the periodic 
"superstructure" are displayed in Fig. 7. It can be seen that the soliton travels slower in 
segments with larger k, and faster in those with smaller k. In the case when the width w 
of each step is relatively small [for instance, w — 1, see panel 7(a)], the soliton evolution 
is clearly non-adiabatic, cf. the above results for the single step (we stress that the full 
modulation period in the present case is essentially large than w, see Fig. 7). Accordingly, 
the soliton gradually develops distortions in its shape and simultaneously emits radiation 
waves. Parallel to distorting itself, the soliton also develops random fluctuations of its 
velocity. 

In the case when the periodic structure is composed of steps with a larger width, for 
instance, w — 10 [see Fig. 7(b)], the soliton also accumulates distortion and generates 
some radiation loss, but, in accordance with the expectation that the evolution must be 
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close to adiabatic in this case, this happens much slower, and the distortion remains mild 
after passing a long distance. With a still larger step width (for instance, w = 15), no 
accumulation of distortion or radiation could be seen at all (the latter case is not shown 
here, as one would simply observe an unperturbed soliton). So, large values of w indeed 
provide for the adiabatic character of the soliton's motion through the periodic structure, 
but the smallest size of w which makes it possible in the case of the periodic system is (quite 
naturally) much larger than the minimum width which provided for the adiabatic passage 
of the gap soliton through the single step, cf. Fig. 6. 

Lastly, it is relevant to mention that virtually identical results were observed in the 
simulations of the periodic system either with periodic boundary conditions (b.c.) in x 
(which corresponds to the above-mentioned case of an apodized BG written in a fiber loop) 
or just in a long domain with absorbing b.c. Thus, the same soliton dynamics is expected 
in both above-mentioned physical realizations of the periodic "Bragg superstructure" . 

VI. CONCLUSION 

In this work, we have studied in detail motion of gap solitons in two models of apodized 
fiber Bragg gratings, including, respectively, a single step of the local reflectivity, or a periodic 
"Bragg superstructure" consisting of alternating steps with opposite signs. Both structures 
offer a potential for the design of various optical elements employing the gap (Bragg) solitons. 
The most important implication of the considered problem is a possibility to halt the soliton 
by passing it through a step with an increasing reflectivity, and eventually to capture a pulse 
of standing light. 

We have developed an analytical approach, which assumes adiabatic evolution of the 
soliton, and is based on the balance equation for its momentum and conservation of the 
energy. The result of the analysis can be obtained in a fully explicit form in the "nonrela- 
tivistic" case, and in an implicit form (as an ordinary differential equation) in the general 
case. Comparison of the predictions produced by the approximation with direct simulations 
shows good accuracy, provided that the width of the inhomogeneity is essentially larger than 
the soliton's intrinsic width. 

Results of systematic direct simulations of the soliton's motion were summarized, showing 
that the soliton is either transmitted across the step or bounces from it, provided that the 
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step is not too narrow; in particular, it is possible to halt the soliton for a very long time. 
If the step is narrow (so that the interaction of the incident soliton with it is no longer 
adiabatic), the soliton may split into two pulses, transmitted and reflected ones (the splitting 
of a gap soliton in a chirped fiber grating has already been observed in the experiment 
We have studied in detail dependences of the outcome of the interaction on the height and 
width of the step, as well as on the initial parameters (velocity and effective mass) of the 
soliton. In particular, a general conclusion is that the soliton is a more cohesive object, 
being less prone to the splitting, if its energy (effective mass) is larger. 

Moving across the periodic structure, the soliton accumulates distortion and radiation 
loss if the structure is composed of narrow steps. The perturbations accumulate much slower 
if the steps are wider, and in the system composed of sufficiently wide steps the soliton can 
move without any loss or irreversible deformation. 
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Figures 



c 




FIG. 1: Borders between the bounce of the gap soliton, initially moving at the velocity cq, from 
the local-reflectivity step with the height of Ak [see Eq. (I14j) ]. and its passage over the step, are 
shown in the parameter plane (Ak, Co), as found from direct simulations, and from the analytical 
approximations - the full ( "relativistic" ) one, based on Eqs. © and and its simplified 

( "nonrelativistic" ) version, which is based on Eq. ()12|) and assumes Cq <C 1 and 9 = const. The 
results are presented for the fixed step's width, w = 20, and fixed mass parameter of the impinging 
soliton, 9 = tt/2. Both Ak and cq are in normalized units. 
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FIG. 2: Different typical outcomes of the collision of the gap soliton with the apodization step 
that has the hight Ak = 0.5 and width w = 6. The upper panels show the waveforms |it(ar)| and 
\v(x)\ (solid and dashed lines) at the end of the simulation, together with the apodization profile, 
n(x). The lower panels depict the evolution of the field in terms of level contours. (All the 

variables plotted are in normalized units.) (a) Deceleration of the soliton with very little radiation 
loss. The initial parameters of the soliton are cq = 0.75 and 6q = 0.57T. 
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FIG. 2: (b) The soliton bouncing back with almost no radiation loss. Here, Co = 0.70 and 
9 = 0.5vr. 
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FIG. 2: (c) An example of the soliton being halted by the step and staying immobile for a long 
time. In this case, cq = 0.7403 and 6q = 0.5n. 
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FIG. 2: (d) Split of the soliton into two pulses, transmitted and reflected ones. In the course 
of subsequent evolution, both pulses gradually decay into radiation. In this case, Co = 0.8 and 
9 = 0.3tt. 



|u|, |v| 

0.9 




FIG. 3: The same as in Fig. 2(d), in the case of a narrower and taller apodization step, with w = 4 
and Ak = 0.7, for the impinging soliton whose mass parameter exceeds the instability threshold, 
9 = 0.6-7T. All the variables plotted are in normalized units. 
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FIG. 4: The share of the initial soliton's energy, which is reflected back by the apodization step 
of the width w = 6, versus the initial velocity c, for different values of the initial mass parameter 
9, and different values of the step's height Ak. Three panels of the figure show the dependence 
e(c) in different regions of c and on different scales, in order to highlight regions where nontrivial 
changes occur. The steep drop of e from I to for larger 9 (such as 9 = 0.4 and 0.5) implies the 
transition from the bounce of the soliton to the transmission without splitting, (e is dimensionless, 
and c is in normalized unit.) 
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FIG. 5: The final velocity c/ of the transmitted pulse and the velocity c& of the reflected one 
(solid lines) are shown, together with the backscattered-energy share e, versus the height of the 
apodization step Ak, for 9 = 0.37T and = 0.57T. The splitting of the incident soliton takes place in 
the interval where both solid curves are present. In this figure, the step's width and initial velocity 
are fixed, w = 6 and cq = 0.8. (e is dimensionless, and cj, q,, Ak are in normalized units.) 
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FIG. 6: The final soliton's velocity Cf (solid line) and the fraction of the backscattered energy, 
e (dashed line), versus the apodized-layer's width w, in the case when the height of the step is 
Ak = 0.5. In this figure, the initial velocity and mass parameter of the soliton are fixed to be 
cq = 0.75 and 6q = 0.5ir. (e is dimensionless, and c/, w are in normalized units.) 
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FIG. 7: Motion of the soliton with the initial velocity Co = 0.9 and initial mass parameter 6q = 0.57T 
through a periodic structure consisting of alternating steps with the local Bragg reflectivity k 
varying between 0.8 and 1. Bottom portions in each panel show the profile of k(x) (dotted line) 
and evolution of the soliton's velocity (continuous line). The top portions show the shape of the 
soliton (|«(x)| and |v(x)|) at the end of the long simulation. (All the variables plotted are in 
normalized units.) (a): The width of the individual step is w = 1 (a non-adiabatic case). 
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FIG. 7: (b): The width of each step is w = 10 (a nearly adiabatic case). 
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